Lattice approach to high-energy hadron-hadron scattering 
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We discuss the non perturbative approach to the problem of high-energy hadron-hadron (dipole-dipole) scat- 
tering at low momentum transfer by means of numerical simulations in Lattice Gauge Theory. 



1. Introduction 

The prediction from first principles of total 
cross sections at high energy is one of the old- 
est open problems of hadronic physics. Present- 
day experimental data are well described by a 
universal pomeron-like power-law behaviour (see, 
for example, Ref. [T] and references therein), 
a tot \ s ) ~ ( s / s o) ep i where the so-called soft 

s — >oo 

pomeron intercept is ep ~ 0.08, but this is forbid- 
den as a true asymptotic behaviour by the well- 
known Froissart-Lukaszuk-Martin theorem [2]. 
As we believe QCD to be the fundamental the- 
ory of strong interactions, it should predict the 
correct asymptotic behaviour; nevertheless, a sat- 
isfactory explanation is still lacking. 

The problem of total cross sections is part of 
the more general problem of high-energy scatter- 
ing at low transferred momentum, the so-called 
soft high-energy scattering. As soft high-energy 
processes possess two different energy scales, the 
total center-of-mass energy squared s and the 
transferred momentum squared t, smaller than 
the typical energy scale of strong interactions 
(1*1 ^5 1 GeV 2 <C s), we cannot fully rely on 
perturbation theory. A genuine non perturba- 
tive approach in the framework of QCD has been 
proposed by Nachtmann in [3] , and further devel- 
oped in |4l5l6l7l8j : using a functional integral ap- 
proach, high-energy hadron-hadron elastic scat- 
tering amplitudes are shown to be governed by 
the correlation function of certain Wilson loops 
defined in Minkowski space. Moreover, as it has 
been shown in |9|10j . such a correlation func- 
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tion can be reconstructed by analytic continua- 
tion from its Euclidean counterpart, i.e., the cor- 
relation function of two Euclidean Wilson loops, 
that can be calculated using the non perturbative 
methods of Euclidean Field Theory. 

In [llj we have investigated this problem by 
means of numerical simulations in Lattice Gauge 
Theory (LGT). Although we cannot obtain an an- 
alytic expression in this way, nevertheless this is a 
first-principle approach that provides (inside the 
errors) the true QCD expectation for the relevant 
correlation function. In this contribution, after a 
quick survey of the non perturbative approach to 
soft high-energy scattering in the case of meson- 
meson elastic scattering, we will present our nu- 
merical approach based on LGT, and we will show 
how the numerical results can be compared to the 
existing analytic models. 

2. Meson meson elastic scattering ampli- 
tudes and Wilson loop correlators 

We sketch here the non perturbative approach 
to soft high-energy scattering; see [Tl] for a more 
detailed presentation. As it has been shown 
in [4|5|6j , the elastic meson-meson scattering am- 
plitude can be reconstructed from the scattering 
amplitude of two qq colour dipoles, after averag- 
ing over the transverse sizes and the longitudinal 
momentum fractions of the dipoles. The central 
quantity in this approach is a certain (properly 
normalised) correlation function (in the sense of 
the QCD functional integral) of two Wilson loops 
in the fundamental representation, defined in 
Minkowski space-time, running along the paths 
made up of the quark and antiquark straight-line 
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classical trajectories and closed at proper times 
±T by straight-line paths in the transverse plane. 

In [9ll0j (see also |12I13| ) it has been shown 
that, under certain analyticity hypotheses, this 
correlation function can be reconstructed from 
the Euclidean correlation function of two Eu- 
clidean Wilson loops, Wi and W2, 



Qe(0;T; 



1,2) = W 



(T)^(T) 



-1 (1) 



(here "1" and "2" stand respectively for Ri±,fi 
and i?2±j /2); the two loops run along two rectan- 
gular paths C\ and C 2 , made up of the "Euclidean 
trajectories" of the partons, 
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where Pipje = ^ ([— ] sin |, Ox, cos - 2 

(o,i? 4± ,o), /? = l-/,, n = (i = 1,2), 

ze = (0, 2_l,0), with and the transverse 
sizes and longitudinal momentum fractions of the 
two dipoles, and z± the impact-parameter dis- 
tance between the two loops in the transverse 
plane. The paths are closed at proper times ±T 
by straight-line paths in the transverse plane; 
here T acts as an IR cutoff which has to be re- 
moved in the end. As the elastic scattering am- 
plitude of two meson states is expected to be an 
IR-finite physical quantity [H], we expect the 
limit T — > 00 to be finite, so that we can define 
C E (0;,gj_;l,2) = G E (0;T -> oo;z±; 1,2). Note 
that is a real function, as can be shown mak- 
ing use of the charge-conjugation invariance of 
the functional integral. 

Finally, the meson-meson scattering ampli- 
tudes can be written as 



M {hh) {s,t;l,2) 



-i 2s 



x J dl |Vi(l)| 2 / d2 |V> 2 (2)| 2 

x C M {0 -> -i\og(s/m 2 ) ;zj.; 1,2);(3) 

here q± is the (transverse) transferred momentum 
(t = — \q±\ 2 ), and ipi an d ^2 are the wave func- 
tions which describe the two interacting mesons. 



Total cross sections are then recovered via the 
optical theorem. 

In the following we will set for simplicity /1 = 
/ 2 = 1/2, which is known to be a good approxi- 
mation for hadron-hadron interactions 11171 . 



3. Wilson loop correlators on the lattice 

The gauge-invariant Wilson-loop correlation 
function Qe is a natural candidate for a lattice 
computation, but some care has to be taken due 
to the explicit breaking of 0(4) invariance on a 
lattice. As straight lines on a lattice can be ei- 
ther parallel or orthogonal, we are forced to use 
off-axis Wilson loops to cover a significantly large 
set of angles. To stay as close as possible to the 
continuum case, the loop sides are evaluated on 
the lattice paths that minimise the distance from 
the true, continuum paths: this can be easily ac- 
complished making use of the well-known Bresen- 
ham algorithm |15j to find the required "minimal 
distance paths'" corresponding to the__sides of the 
loops. The relevant Wilson loops Wl(Z|| ; r±; n) 
are then characterised by the position n of their 
center and by two two-dimensional vectors lu and 
r_i_ , corresponding respectively to the longitudinal 
and transverse sides of the loop. 

On the lattice we then define the correlator 



GL(h\\,k\\;d±',f f i±,f f 2±) 

" (W L (h || ; rU 5 d)) (W L (f 2 „ ; f 2X ; 0)) 



where d — (0, d±,0), d± = (<i 2 , (Z3), and moreover 
CL(h\\,hfd±;r 1 ±,r 2 ±) 

ee lim ^iGi||,/ 2 ||;d_L;fix,f 2 j_),(5) 

Li,L,2 — >oo 

where Li = \lq | arc defined to be the lengths of 
the longitudinal sides of the loops in lattice units, 
and ij|| ee Lili\\. In the continuum limit, where 
0(4) invariance is restored, we expect 

CL(h\\,l2\\',dx;fi±,f2±) 

~ C E {0; ad±;ar 1± , 1/2, ar 2 ±, 1/2), (6) 

where lu ■ Z 2 || = cos# defines the relative angle 8, 
and a is the lattice spacing. 
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Figure 1. Lattice data plotted against for vari- 
ous lengths of the loops. 

4. Numerical results 

As already pointed out in the Introduction, nu- 
merical simulations cannot provide the analytic 
expression for the relevant correlation function, 
but nevertheless, as these simulations are first- 
principles calculations, they provide the "correct" 
(inside the errors) prediction of QCD. Approxi- 
mate analytic calculations have then to be com- 
pared with the lattice data, in order to test the 
goodness of the approximations involved. In par- 
ticular, we are interested in the dependence on 
the relative angle 9, as it encodes the energy de- 
pendence of the scattering amplitudes, which is 
recovered after the proper analytic continuation. 

In Fig. Q] we show, as an example, the lattice 
data for Gl hi the case of parallel transverse sides 
with |r*ij_ | = |r2j_| = 1 at d = 0, plotted against 
the angle 9 for various lengths of the loops. These 
data are obtained using Wilson action for SU(3) 
pure-gauge (quenched) theory, on a 16 4 lattice 
at f3 = 6.0. The data are quite stable against 
variations of the lengths, so that we can take the 
largest-length data as a reasonable approxima- 
tion of Cl- 

In Figs. [21 and [3] we compare Cl with the pre- 
diction of various models (the loop configuration 
is the same as in Fig. [1]). While the Stochas- 
tic Vacuum Model (SVM) [TB] provides a fully 
quantitative prediction, that can be directly com- 
pared with the data, the Instanton Liquid Model 
(ILM) [17] and the AdS/CFT correspondence [18] 
give only the qualitative dependence on the an- 
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Figure 2. Comparison of lattice data with the 
SVM prediction (solid) and with a best-fit with 
the SVM functional form (dotted). 
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Figure 3. Comparison of lattice data with best- 
fits with the lowest-order perturbative (solid), 
the ILM (dotted) and the AdS/CFT (dashed) ex- 
pressions. 

gle 9, so that a comparison can be made by 
trying to fit the data with the given functional 
form. In Fig. [2] we show the SVM prediction, to- 
gether with a best-fit with the SVM functional 
form; in Fig. [3] we show the best-fits with the ex- 
pressions obtained in perturbation theory to low- 
est order [19116110] . in the ILM and using the 
AdS/CFT correspondence. A detailed discussion 
of the results is given in [11) ; here we simply note 
that the agreement of the numerical data with the 
various models is not fully satisfactory, and fur- 
ther investigations have to be made, both on the 
numerical and on the analytical side. We want 
also to remark that while perturbative effects 



Figure 4. Dependence of Ce on the distance for 
9 = 45° and 9 = 90° in the case ri± = f 2 j_, 
|ri_i_| = 1, ri± || dj_ (logarithmic scale). 




Figure 5. "Antisymmetric" part of lattice data, 
and corresponding SVM prediction. 

seem to be dominant at short distances between 
the loops, non perturbative effects are already rel- 
evant at distances of about 0.2 fm; however, as the 
correlation function is rapidly decreasing with the 
distance d = \dj_ \ between the centers of the loops 
(see Fig. |4]), a detailed study at large distances is 
difficult, and requires the use of noise reduction 
techniques. 

Lattice data show also the presence of odd- 
eron contributions to dipole-dipole scattering. 
Indeed, as explained in [13111] . making use of 
the crossing-symmetry relations for loops p~2] one 
can show that the crossing-odd component of the 
dipole-dipole scattering amplitudes is related via 
the usual analytic continuation to the antisym- 
metric (with respect to 7r/2) component of 
Ce- this quantity is shown in Fig. [5j together 



with the coresponding SVM prediction (the loop 
configuration is the same as in Fig.Q]). 
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